QED Effective Actions in Inhomogeneous 

Backgrounds: 
Summing the Derivative Expansion 



Gerald V. Dunne 

O 

', Department of Physics, University of Connecticut, Storrs, CT 06269-3046, USA 

(N 

> 
O 

. Abstract. The QED effective action encodes nonlinear interactions due to quantum 

vacuum polarization effects. While much is known for the special case of electrons 
in a constant electromagnetic field (the Euler-Heisenberg case), much less is known 
for inhomogeneous backgrounds. Such backgrounds are more relevant to experimental 
^ . situations. One way to treat inhomogeneous backgrounds is the " derivative expansion" , 

' in which one formally expands around the soluble constant-field case. In this talk I 

. use some recent exactly soluble inhomogeneous backgrounds to perform precision tests 

on the derivative expansion, to learn in what sense it converges or diverges. A closely 
related question is to find the exponential correction to Schwinger's pair-production 
' formula for a constant electric field, when the electric background is inhomogeneous. 



This talk is concerned with the one-loop QED effective action [1,2]: 
r^/- S[A] = -^\ogdet[p^ + 



I Here = dfj_ + ieA^ is the covariant derivative, and so S[A] is a functional of the 

^ ' classical background field A^lx). The effective action is the generating functional 

for one-fermion-loop Green's functions, which describe the nonlinear QED effects 
to this order. Ideally, one would like to know S[A] for any field A^{x), but this is 
not feasible. However, for the special case where the field strength tensor, F^,^ = 
dfj^Ai, — d^Afj^, is uniform, the effective action can be evaluated in closed form [3-5]: 

S = ^{E' - B') + / d'x \{E' - By + 7{E ■ Bf] +... (2) 

Here the first term is the familiar classical Maxwell action, while the next term 
gives the leading quantum correction, which is quartic in the field strengths. The 
fine structure constant a = |^ in these units. All higher terms in the expansion (2) 
are known explicitly [3-5]. It can also be shown that when the background field 
is a constant electric field of strength E, the effective action has an exponentially 
small imaginary part 
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This imaginary part of S gives the pair-production rate for electron-positron pro- 
duction from vacuum in a strong uniform electric field [3,5,6]. Note also that it is 
manifestly nonperturbative in the coupling constant e, in contrast to the perturba- 
tive series expansion in (2). 

For a general inhomogeneous background field we need some sort of approxima- 
tion to determine the effective action. One such approximation is the derivative 
expansion [7-10], in which one assumes that F^,^ is "slowly varying", so that one 
can make the following formal expansion 

S = [F] + 5(2) [F, {dFf] + ... (4) 

where S''^^"-' contains 2n derivatives of Fn^,. Clearly, this is a rather formal and ill- 
defined expansion. In this talk, I ask two specific questions about this expansion: 

1. Does the derivative expansion converge? 

2. How is the nonperturbative imaginary piece in (3) modified? 
These two questions are in fact intimately related. 

First, let us study the convergence properties of the constant-field Euler- 
Heisenberg effective action (2). When the uniform background is purely magnetic, 
the entire series (we ignore the Maxwell term) in (2) is 



7r2 VmV ^0 (2n + 4)(2n + 3)(2n + 2) VmV 

This is a " low energy" effective action in the sense that we assume that the char- 
acteristic cyclotron energy of the magnetic background, is much less than the 
electron rest energy, m. The expansion coefficients in (5) involve Bernoulli numbers 
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[11] alternate in sign and grow factorially in magnitude. In fact, for large n, 
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Thus, the Euler-Heisenberg effective action is itself a divergent series! This is not 
a bad thing; perturbation theory is typically divergent. In many examples (the 
anharmonic oscillator, the Stark effect, 0^, 0^, QED, ...) perturbation theory is 
known to be divergent [12]. Typically, the expansion coefficients diverge as 

|a„| -/3"r(7n + 5)[l + ...] (7) 



at large orders n, and usually the parameter 7 = 1 or 2. 

Even though a series may be divergent, we can still extract important and mean- 
ingful information from it. One such approach is known as Borel summation [13,14]. 



To motivate this approach, consider the following expression for the (clearly diver- 
gent) alternating series 
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Read from right to left, (8) says that the asymptotic expansion for small g of the 
integral on the RHS (which converges for all g positive) yields the series on the 
LHS. Read from left to right, (8) involves a formal interchange of summation and 
integration (a step whose validity must be examined carefully in a given case). We 
call the integral on the RHS the Borel sum of the divergent series on the LHS. If 
we were to try this for a non- alternating series, we might write analogously 

f{-9) = Y.n\9--- I ds-— (9) 

n=0 9 ■'0 i. b 

This is problematic, as there is a pole at s = 1 on the integration contour. This pole 
must be resolved, and in doing so, an imaginary part arises. Roughly speaking, this 
imaginary part is of the form /m/(—5f) = ^ exp[— -^J, which is non-perturbative in ^f. 
Similar arguments apply if the expansion coefficients grow like |a„| ~ /5" r(7 n + 5), 
rather than just n\ as in the above example. Then f{g) = Y,n^n9^ behaves like 





Imf{-g)^-\ — \ exp[- — ] (11) 



I caution that these Borel representations are somewhat formal, as we are assuming 
that there are no additional poles or cuts in the complex g plane to prevent the 
straightforward analytic continuation on which these expressions are based [13,14]. 
In this talk our attitude is exploratory; we shall use these formulae cautiously, with 
cross-checks, but defer questions such as uniqueness to more rigorous studies. 

It is an instructive exercise to apply this Borel technology to the alternating 
divergent Euler-Heisenberg effective action (5) [15-17]. Combining the asymptotic 
form (6) of the coefficients with the Borel formula (10) we obtain 

e^B'^ r°° ds 1 s\ , ni?s, 
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where we have used: coths - 1/s = 2s/7r^E'^^i l/(/^^ + sV^r^), and E^i = 
7r^/6. But this is precisely Schwinger's proper-time integral representation [5,3,4] of 
the effective action, which we see can be viewed as the Borel sum of the (divergent) 
Euler-Heisenberg perturbation series (5). If the constant background field is electric 



instead of magnetic, the only change perturbatively is to replace B'^ in (5) by —E'^ 
(because the Lorentz invariant combination is E'^ — B'^). Thus the alternating 
divergent series (5) becomes a non- alternating divergent scries. Applying the Borcl 
formula (11), together with the asymptotic growth rate in (6), we find an imaginary 
part in complete agreement with Schwinger's result (3). So, in this case, our use of 
the Borel formulae (10-11) is consistent. Further, we see that the Euler-Heisenberg 
series (5) had to he divergent. If it were not, then there would be no essential 
difference between the magnetic and electric cases, and we would not find any 
imaginary part for the electric case; so we would miss the genuine physical effect 
of vacuum instability. This is somewhat reminiscent of Dyson's physical argument 
[18] that QED perturbation theory should be non-analytic as a series in a = 
since a negative is unstable. 

Unfortunately, the pair production rate derived from (3) is extremely small. 
Indeed, the critical electric field, E^. = ~ 10^^ V/ cm, where the exponent is 
of order 1, is well beyond laboratory static fields. However, using intense laser 
fields, it has recently become possible to probe this critical regime [19]. These 
involve inhomogeneous fields, and so it becomes important to ask how the Euler- 
Heisenberg analysis is modified by a field inhomogeneity. Thus we return to our 
two questions concerning the convergence properties of the derivative expansion (4) 
and how this might modify the imaginary part (3). 

A first guess is that one should look at the high orders of the derivative expansion 
to see if it is diverging. Unfortunately this is impossible for two fundamental 
reasons. First, the derivative expansion is not actually a series. This is because as 
one goes to higher orders, there is a rapid proliferation of new terms that do not 
have a counterpart at previous orders [8-10]. Second, it is extremely difficult to 
calculate any of these terms at very high order, so we could not expect anyway to 
observe any asymptotic behaviour of the expansion coefficients. 

However, we can avoid both these problems at once by considering some special 
exactly solvable cases [15]. In these cases, the background inhomogeneity can be 
characterized by a single scale parameter, so that the derivative expansion becomes 
a true series, in inverse powers of this scale parameter. And since these cases are 
soluble, we have access to all orders of the derivative expansion, so that we can 
probe the high order behaviour. So, while we sacrifice generality by concentrating 
on specific backgrounds, we gain the ability to perform precise analytic tests of the 
derivative expansion. 

To be specific, we will use the fact that for the following backgrounds, the effective 
action has been computed in closed-form [20,21]. 

B{x) ^ Bsech^i^) or E{t) ^ E sech^i-) (13) 

A T 

It has of course been known for many years that the corresponding Dirac equations 
are soluble [22-24]. But it is still non-trivial to perform exphcitly the necessary 

traces so as to express the effective action (1) as a series or as an integral repre- 
sentation, in terms of a single integral, just as in the Euler-Heisenberg case (5) or 



(12). For the inhomogeneous magnetic background, B[x) = 
the exact effective action is [20] 
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Several comments are in order concerning this result. First, note that the se- 
ries expansion in (14) is a double sum, with derivative expansion parameter, 
and perturbative expansion parameter, Second, the expansion coefficients are 
known to all orders, and are relatively simple numbers, just involving the Bernoulli 
numbers and factorial factors. Third, it has been checked in [20] that the first 
few terms of this derivative expansion are in agreement with explicit field theoretic 
calculations, specialized to this particular background. Fourth, a simple integral 
representation for (14) can be found in [20]. 

Given the explicit series representation in (14), we can check that the series 
is divergent, but Borel summablc (in the sense of our earlier discussion), in the 
magnetic case. This can be done in several ways. One can either fix the order k 
of the perturbative expansion in (14) and show that the remaining sum is Borel 
summable, or one can fix the order j of the derivative expansion in (14) and show 
that the remaining sum is Borel summable. Or, one can sum explicitly the k sum, 
for each j, as an integral of a hypergeometric function, and show that for various 



values of the remaining derivative expansion is divergent but Borel summable. 
These arguments do not prove rigorously that the double series is Borel summable, 
but give a strong numerical indication that this is the case. 

The case of the inhomogeneous electric field, E(t) = £^sech^(^), in (13) can also 
be solved explicitly [21]. A short-cut to the answer is to note that we can simply 



make the replacements, B'^ 



-E\ and A' 



-T , in the magnetic case result 



(14). In particular this has the consequence that the alternating divergent series 
of the magnetic case becomes a non- alternating divergent series, just as was found 
in the Euler-Heisenberg constant-field case. For example, fixing the order j of the 
derivative expansion, the expansion coefficients behave for large k (with j fixed) as 



^ ^y^k V{2k + j)V{2k + J + 2)g2fc+2,+2 ^ , r(2fc + 3j-|) 
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Note that these coefficients are non-alternating and grow factorially with 2k, as in 
the form of (7). Applying the Borel dispersion formula (11) gives 
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Remarkably, this form can be resummed in j, yielding a leading exponential form 



We recognize the first term in the exponent as the famihar Schwinger exponent from 
(3), and so the second term may be viewed as the leading exponential correction 
to the constant-field answer (3). This is what wc set out to find, and we see 
that it arose through the divergence of the derivative expansion. I stress that this 
exponential correction is not accessible from low orders of the derivative expansion, 
such as those studied in [8-10]. 

But the situation is even more interesting than this result (17) suggests. For 
example, we could instead have considered doing the Borel resummation for the j 
summations, at each fixed k. Then for large j, the coefficients behave as [15] 
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which once again are non-alternating and factorially growing. Applying the Borel 
dispersion formula (11) gives 

,,)^^(2^eEr^ (19) 
4^3^5/2 (2A;)! ^ ' 

Once again we see that this leading form can be resummed, leading to 

777,3/2 r ^ eET\' 

I-^-^;^-p[-2.mr(l-^)J (20) 

But this leading exponential form of the imaginary part is different from that 
obtained in (17), and moreover, it is different from the constant-field case (3). So, 
what is going on? The answer is that there are two competing leading exponential 
behaviours buried in the double sum (14), and the question of which one dominates 
depends crucially on the relative magnitudes of the two expansion parameters. 
These two expansion parameters are the derivative expansion parameter, and 
the perturbative expansion parameter, Another important parameter is their 
ratio, since this sets the scale of the corresponding gauge field: 

A{t) _ eETtanh(t/T) eEr _ eE/w? 

(21) 



m m m l/(mT) 



Thus, we can define a "non-perturbative" regime, in which ^ > 1. Then 

2 2 

rriT 3> so that the dominant exponential factor is cxp[— ^] ^ exp[— 27rmr]. 
In this regime, the leading imaginary contribution to the effective action is given 
by the expression (17), and we note that it is indeed non-perturbative in form, 
and the correction in the exponent is in terms of the small parameter <^ 1. 
Alternatively, we can define a "perturbative" regime, in which ^ <^ 1. Then 
mr <^ so that the dominant exponential factor is exp[— 27rmr] ^ exp[— ^]. In 
this regime, the leading imaginary contribution to the effective action is given by 



the expression (20), and we note that this is in fact perturbative in nature (despite 
its exponential form). 

To understand these two regimes more carefully, we can use the WKB approach 
developed by Brezin and Itzykson [25] , who studied the case of a sinusoidal electric 
field: E{t) = Ecos{u!t). This case is not exactly soluble, but a WKB expression 
for the imaginary part of the effective action is: 

ImS ~ / d^k exp[-7rQ] (22) 



where fl = ^ j^^'^m? + k\ + [k^ — eAz{f)Y. If wc apply this same WKB anal- 
ysis to the (exactly soluble) case E{t) = Escch^^t/r) in (13) we find that 
n = T{^m? + kl + {cEt + k^y + + k'i + {eEr - k^f - 2eEr). Moreover, 
this system is known to be WKB-exact [26,24]. Then when the momentum inte- 
grals in (22) are done, we obtain precisely the leading results (17) or (20), depending 
on whether we are in the "non-perturbative" ^ » 1, or "perturbative" ^ <^ 1 
regime [15]. This serves as a useful cross-check of our Borel analysis. 

To conclude, there exist two inhomogeneous backgrounds (one magnetic and one 
electric) for which the effective action (1) can be evaluated exactly, in closed- form, 
thereby generalizing the well-known Euler-Heisenberg constant-field case. In this 
talk I have used these cases to probe the convergence properties of the derivative 
expansion. We learn that the derivative expansion is divergent. In the magnetic 
case it appears to be Borel summable, and in the electric case the direct application 
of Borel dispersion relations leads to imaginary parts that are in complete agreement 
with an independent WKB analysis. Indeed, we can argue a la Dyson that the 
derivative expansion must be divergent, otherwise there would be no exponential 
correction to the imaginary part of S for an inhomogeneous electric field. There 
are two competing leading exponentials, and in the non-perturbative regime, where 
3> 1, we find an exponential correction (17) to Schwinger's formula (3). In the 
perturbative regime, where <^ 1, a different exponential factor dominates. 



REFERENCES 

1. W. Dittrich and M. Reuter, Effective Lagrangians in Quantum Electrodynamics, 
(Springer, Berlin, 1985). 

2. W. Greiner and J. Reinhardt, Quantum Electrodynamics, (Springer, Berlin, 1992). 

3. W. Heisenberg and H. Euler, "Folgerungen aus der Diracschen Theorie des 
Positrons", Z. Phys. 98, 714 (1936). 

4. V. Weisskopf, "Uber die Elektrodynamik des Vakuums auf Grund dcr Quantcn- 
theorie des Elektrons", Kong. Dans. Vid. Selsk. Math-fys. Medd. XIV No. 6 (1936), 
reprinted in Quantum Electrodynamics, J. Schwinger (Ed.) (Dover, New York, 1958). 

5. J. Schwinger, "On Gauge Invariance and Vacuum Polarization", Phys. Rev. 82, 664 
(1951). 



6. F. Sauter,"Uber das Verhalten eines Elektrons im homogenen elektrischen Feld nach 
der relativistischen Theorie Diracs", Z. Phys. 69, 742 (1931). 

7. I. J. R. Aitchison and C. M. Eraser, "Derivative Expansions of Eermion Determi- 
nants: Anomaly Induced Vertices, Goldstone-Wilczek Currents and Skyrme Terms" , 
Phys. Rev. D 31, 2605 (1985). 

8. D. Cangemi, E. D'Hoker and G. Dunne, "Derivative Expansion of the Effective 
Action and Vacuum Instabihty for QED in 2 + 1 Dimensions", Phys. Rev. D 51, 
R2513 (1995). 

9. V.P.Gusynin and I.A.Shovkovy, "Derivative Expansion for the One-Loop Effective 
Lagrangian in QED", Can. J. Phys. 74, 282 (1996). 

10. D. FHegner, P. Habcrl, M.G. Schmidt, and C. Schubert, "The Higher Derivative 
Expansion of the Effective Action by the String Inspired Method (Part 2)", Ann. 
Phys. (NY) 264, 51 (1998). 

11. The BernouUi numbers B2n are defined by x/{e^ — 1) = J2'^=q BnX"'/n\. 

12. For an excellent review, see: J. C. Le Guillou and J. Zinn-Justin (Eds.), Large-Order 
Behaviour of Perturbation Theory, (North Holland, Amsterdam, 1990). 

13. G. Hardy, Divergent Series (Oxford Univ. Press, 1949). 

14. C. M. Bender and S. A. Orszag, Advanced Mathematical Methods for Scientists and 
Engineers (McGraw-Hill, New York, 1978). 

15. G. Dunne and T. Hall, "Borel summation of the derivative expansion and effective 
actions", Phys. Rev. D 60, 065002 (1999). 

16. A. Zhitnitsky, "Is an Effective Lagrangian a Convergent Series?" , Phys. Rev. D 54, 
5148 (1996). 

17. G. Dunne and C. Schubert, "Two-loop Euler-Heisenberg QED pair-production rate", 
Nucl. Phys. B 564, 591 (2000). 

18. F. J. Dyson, "Divergence of Perturbation Theory in Quantum Electrodynamics", 
Phys. Rev. 85, 631 (1952). 

19. D. Burke et al, "Positron Production in Multiphoton Light-by-Light Scattering", 
Phys. Rev. Lett. 79, 1626 (1997); A. Melissinos, "The Spontaneous Breakdown of 
the Vacuum", hep-ph/9805507. 

20. G. Dunne, T. Hall, "An Exact QED3+1 Effective Action", Phys. Lett. B 419, 322 
(1998); D. Cangemi, E. D'Hoker, G.Dunne, "Effective Energy for QED2+1 with Semi- 
localized Static Magnetic Field: A Solvable Model", Phys. Rev. D 52, R3163 (1995). 

21. G. Dunne and T. Hall, "QED effective action in time dependent electric back- 
grounds", Phys. Rev. D 58, 105022 (1998). 

22. F. Sauter, "Zum 'Kleinschen Paradoxon"', Z. Phys. 73, 547 (1932); this includes a 
related solvable Dirac equation for a spatially inhomogeneous electric background. 

23. See, e.g., N. Narozhnyi and A. Nikishov, "The Simplest Processes in a Pair-producing 
Electric Field", Sov. J. Nucl. Phys. 11, 596 (1970), and refs. therein. 

24. A. B. Balantekin, J. E. Soger and S. H. Frickc, "Dynamical Effects in Pair Production 
by Electric Fields", Int. J. Mod. Phys. A 6, 695 (1991). 

25. E. Brezin and C. Itzykson, "Pair production in Vacuum by an Alternating Field", 
Phys. Rev. D 2, 1191 (1970); E. Brezin, These (Paris, 1970). 

26. A. Comtet, A. Bandrauk and D. Campbell, "Exactness of Semiclassical Bound State 
Energies for Supersymmetric Quantum Mechanics", Phys. Lett. 150 B, 159 (1985). 



